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1. Introduction 

Radiative kaon decays offer a unique possibility to study Chiral Perturbation Theory (ChPT) in 
detail. In particular, direct photon emissions as in K"^ 71"'' 71° 7 decays or in decays with vanishing 
^{p^) as 71* 77'*' are of theoretical interest. 

The NA48/2 experiment has collected data on charged kaon decays in the years 2003 and 
2004. The kaon beams had a momentum of 60 GeV/c with and K decays being recorded 
simultaneously, to systematic effects in CP violation measurements [|l]]. The data were recorded 
with both a highly efficient 3-track-trigger for decays of charged kaons into three charged particles, 
and a 1 -track-trigger, which required a minimum invariant mass of the neutral decay particles to ex- 
clude the abundant 71* 7r" and — > /x*v^ decays. In total, several billions of reconstructed 
decays were recorded. 

The NA48 detector is described in detail elsewhere [^. The main detector components were 
a magnetic spectrometer, consisting of two sets of two drift chambers before and after a dipole 
magnet, providing a momentum resolution of about 1 .4% for 20 GeV/c charged tracks, and a liquid- 
krypton electromagnetic calorimeter (LKr) with an energy resolution of about 1% for 20 GeV 
photons and electrons. 



2. :;r*:;r^ 7 Decays 

The total amplitude of the 71* 7r'' 7 decay is the sum of two terms: inner bremsstrahlung 

(IB), with the photon being emitted from the outgoing charged pion, and direct emission (DE), 
where the photon is emitted from the weak vertex. The IB component can be predicted from QED 
corrections to A'* 71*71° in a straight-forward way [||, ^. For the DE term, several studies within 
the framework of Chiral Perturbation Theory (ChPT) exist [|, |, 0, |, ^. At (ff{p^) ChPT, direct 
photon emission can occur through both electric (X^) and magnetic (X^) dipole transitions. The 
magnetic part is the sum of a reducible amplitude, that can be calculated using the Wess-Zumino- 



Witten functional |10, 11 1, and a direct amplitude, which size is expected to be small. For the 



electric transition no definite prediction exists. 

The total decay rate as a function of the kinematic variable = {p^ ■ Py){pK • Py) / {^K^n) 
is given by 

ar* dvt 



IB 



[l+2cos(±0 + 5/-52)|X£|w2 + m>i(|X£p + |XMp)W^]. (2.1) 

In addition to the IB and DE contributions, the decay rate contains also the interference (INT) 
between IB and DE, which, apart of the strong 7171 re-scattering phases 5/ and 5q , depends only on 
Xe and a possible CP violating phase 0. By measuring the INT term it is possible to disentangle 
the electric and magnetic amplitudes and to investigate possible CP violation in — > 7r*7r*'7. 

Previous measurements have been performed by several experiments. The combined DE 
branching fraction, based on the world total of about 30000 /T* 7r*7r*'7 events, is Br(DE) = 
(4.3 ±0.7) X 10^^ 1 12], with the assumption of no interference term, consistent with the only previ- 
ous measurement of Frac(INT) = Br(INT) /Br(IB) = (-0.4± 1.6)% by the E787 experiment [0]. 
All previous measurements were performed in the restricted kinematic region 55 < < 90 MeV 
of the pion kinetic energy in the kaon rest frame. 
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NA48/2 is the first experiment wliicli can use botii and K events. In addition, a strong 
suppression of — > n'^n^n^ events, based on the excellent performance of the LKr calorimeter, 
was implemented. This allowed to extend the kinematic region to < < 80 MeV, with a slightly 
stronger upper cut due to the on-line trigger rejection of K"^ — > tt^tt" events. The remaining back- 
ground, coming mainly from K"^ n^n'^n^, was estimated with Monte Carlo simulated events 
to be less than 1 % of the DE contribution. The probabihty of mis-identifying the odd photon was 
estimated to be less than 10^^. 




Figure 1: Left: Selected K tt tt / candidates. Center: Maximum-likeUhood fit of tfie W distribution 
of the selected ;r*;r''7 candidates. Right: Fit residuals. 



In total, about 1 million of —>■ tt^tt"/ events were reconstructed by NA48/2 (Fig. |I] (left)). 
The extraction of the IB, DE, and INT contributions was done with an extended maximum-likelihood 
fit of the Monte Carlo W distributions of the single components to the data distribution. For 
the fit, the gamma energy was required to be above 5 GeV to be insensitive of inefficiencies 
of the LI trigger for small cluster energies. In addition, the kinematic range was restricted to 
0.2 <W < 0.9, leaving about 600000 events for the fit. The fit to the data is shown in Fig. H and 
yielded Frac(DE) = (3.32 ±0.15)% andFrac(INT) = (-2.35 ±0.35)% (for < < 80 MeV). 

As a cross-check, a simple polynomial fit to the data W distribution, divided by the Monte 
Carlo IB distribution was performed (Fig. ^ (left)). Although this method does not fully correctly 
take into account the acceptances, the result was in perfect agreement with the maximum-likelihood 
method. 




Figure 2: Left: Polynomial fit to the ratio of data over IB Monte Carlo. Right: Contour plot for the DE and 
INT terms. The black cross shows the lc7 statistical uncertainties of the projections. 
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Many possible systematic uncertainties were investigated. Most contributing were the descrip- 
tion of the detector acceptance, trigger efficiencies, and the LKr calorimeter energy scale. The final 
result, including also systematic uncertainties, is 

Frac(DE)o<7-,*<80MeV = ( 3.32 ±0.15stat ±0.1 Vt) x IQ-^, (2.2) 
Frac(INT)o<r*<8()Mev = (-2.35 ±0.35stat ±0.39syst) x lO^^, (2.3) 

with a correlation coefficient of —0.93 between both values. Fig. ^ (right) shows the confidence 
regions for the statistical uncertainties. 

From this, the electric and magnetic amplitudes can be extracted to 

Xe = (-24 ± 4stat ± 4,yst) GeV-^ (2.4) 
Xm = (254± llstat ± llsyst) GeV-^ (2.5) 

with the magnetic amplitude being very close to the WZW prediction of about 271 GeV^"* [^, p^ ]. 

For comparison with previous experiments, a fit with the ESJT term set to was performed. 
The result, extrapolated to the kinematic range 55 < T* < 90 MeV, was 

Br(DE)iNT=?^90M^v = (2.32 ± 0.05,tat ± 0.08,yst) x 10-^ (2.6) 

in clear disagreement with the previous measurements. The X'^ of this fit was 51.0/12 (compared 
to 14.3/11 when including the INT term as a free fit parameter), strongly indicating the need of the 
INT term for a proper description of the data. 

Finally, possible direct CP violation in this channel was investigated. CP violation would 
manifest itself in a decay rate asymmetry of with respect to K decays and/or in different W 
distributions for and K , due to a non- vanishing phase ^ in the differential decay rate. A 
possible decay rate asymmetry can be expressed in an asymmetry of the total number of events, 
defined as An = {N+ - RN-) / + RN-), with and the numbers of and K decays, 
and R the ratio of to K in the beam, determined from K"^ n^n'^n^ decays ^ Removing the 
cuts on the W range and the photon energy, thus using the complete data set of more than a million 
decays, NA48/2 found A/v = (0.0± l.Ostat ±0.6syst) x 10"^, corresponding to \An\ < 1.5 x 10"^ at 
a confidence level of 90%. Extraction of the CP violating phase (j) yielded sin(j) = —0.01 ±0.43, 
equivalent to | sin0| < 0.56 at 90% CL. 

Assuming the interference term to be the origin of possible CP violation, a fit to the ratio of 
the W spectra of and K , given by ^ = (l + (a it e)W^ + bW^) , was performed. With 
the parameters a and b from the DE and INT fractions, a single parameter fit obtained Aw = 
e J (INT/IB) = (-0.6 ± 1.0) x 10"\ in good agreement with the value of A 



N- 



3. :;r^ 77 Decays 

The —>■ n^YY decay is of high interest in ChPT, since contributions of &{p^) vanish, thus 
giving high sensitivity to <ff{p^) and &{p^). The differential decay rate of —>■ 71^77 is given as 



dydz 2%3 



z2(|A + B|2+|C|2) + (/_!a(1,4,^)) (|Bp + |D|2) 



(3.1) 



'This assumes negligible CP violation in n^n^n^^, which is consistent with the NA48/2 limit on CP violation 

in the k^k^k^ Dalitz plot JII]. 
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with y = (£■* -E^)/mK and z = m^^/m^. 

At ff{p^), predominantly loop diagrams contribute, leading to a distinct cusp in the invariant 
77 mass at twice the 71+ mass [0]. The amplitude A is known up to a parameter c, which needs 
to measured from experiment. Also at i^{p^), poles and tadpole diagrams contribute to the C 
amplitude [ 15 1 . At ^ {p^ ) , unitarity corrections could alter the branching fraction by 30 — 40% [ ]T^ ] . 

So far analyzed were about 40% of the complete data set. Due to the similarity in topology 
to K"^ 71+71° events, which were trigger suppressed, the signal trigger efficiency was only about 
40%. In total 1164 candidates were reconstructed and passing the selection, corre- 

sponding to about 40 times the previous world sample. The background contribution, mainly from 
71+71*' 7 events, was determined from Monte Carlo simulation to 3.3%. The invariant 71+77 
and 77 mass distributions are shown in Fig. ^, the latter exhibiting the expected cusp at twice the 
pion mass. 





0.32 0.34 0.36 
m^., [GeV/c^] 



Figure 3: Selected — + tt 77 candidates. Left: Invariant tt 77 mass. Right: Invariant 77 mass. 

Obtaining the detector acceptance from a simulation using ff{p^) ChPT with c = 2, a prelim- 
inary, model-dependent branching fraction was obtained: 

Br(/:+ ^ 7r±77),=2,^(/) = (1.07 ±0.04,tat ±0.08syst) x 10"^ (3.2) 

The systematic uncertainty is dominated by the trigger efficiency. A model-independent measure- 
ment and the extraction of the parameter c are in preparation. 



^.K^^n^e^e /Decays 

The decay —>■ 7r+e+e^7 is similar to K"^ — > 71+77, with one of the photons internally 
converting into a pair of electrons. As for K"^ —>■ n^yy, in ff{p^) ChPT the branching fraction and 
the e+e^7 spectrum are determined by a single parameter c. At &{p^), unitarity corrections may 



alter the branching fraction by up to 40% [17]. 

Using the whole NA48/2 data set, 120 signal candidates with only small background conta- 
mination were found (Fig. ^. This is the first observation of this decay. 

As normalization channel the abundant decay K"^ —> n'^n^ with — > e+e^7 was used. The 
branching fraction was computed in bins of nieey, thus being independent of any assumption on the 
meey distribution. Integrating over the single bins in the accessible region gave [ 18 1 

Bt{K^ ^ 7i^e+e-y)„^^^^260MeV/c^- = ( 1. 19 ±0. 1 2,t,t± 0.04,^,0 x 10"^ (4.1) 
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. Data 

I I K*-> 7I*II°II° 

I I K*^ 7i*ii°e*e" 
I I K*-> 7i*e*e" 




0.45 0.5 0.55 

Invariant Ji*e*e"Y mass [GeV/c^] 




0.2 



0.25 0.3 0.35 

Invariant e*e"y mass [GeV/c^] 



Figure 4: Invariant TT+e+e 7 (left) and e+e /mass for the selected /T''^ ^ TT+e+e /candidates. 

A single-parameter fit to the nieey distribution above 260 MeV/c^ gave a value c = 0.90 ± 0.45. 
Using tliis value for c, the total branching ratio was obtained as 



Br{K^ Ti^e^e'r) = (1.29 ±0.13exp ±0.03^) x 10"^ 
where the last uncertainty reflects the model uncertainty for mgey below 260 MeV/c^. 
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